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Abstract 

The compressibility equation of state for a multicomponent fluid of particles 
interacting via an infinitely narrow and deep potential, is considered within 
the mean spherical approximation (MSA). It is shown that for a class of 
models leading to a particular form of the Baxter functions qij(r) containing 
density-independent stickiness coefficient, the compressibility EOS does not 
exist, unlike the one-component case. The reason for this is that a direct 
integration of the compressibility at fixed composition, cannot be carried out 
due to the lack of a reciprocity relation on the second order partial derivatives 
of the pressure with respect to two different densities. This is, in turn, related 
to the inadequacy of the MSA. A way out to this drawback is presented in 
a particular example, leading to a consistent compressibility pressure, and a 
possible generalization of this result is discussed. 



I. INTRODUCTION 



Baxter's 'sticky hard sphere' model 0],Q] (hereafter referred to as SHS1 model) has often 
been employed in studies on colloidal suspensions of adhesive particles. Its potential adds 
to a hard sphere (HS) repulsion an infinitely strong surface adhesion, defined by taking 
an attractive square-well tail with vanishing width and infinitely increasing depth, giving 
a finite non-zero contribution to the second virial coefficient ('sticky limit') The SHS1 
model admits analytical solution if the Ornstein-Zernike (OZ) integral equations of the liquid 
state theory are coupled with the Percus-Yevick (PY) approximation |T|-^[. The resulting 
expression for a fluid with p components requires the knowledge of a set of density-dependent 
parameters {A^}, whose values have to be determined by solving numerically p(p + l)/2 
coupled quadratic equations @. The applicability of the SHS1-PY solution is, therefore, 
limited to systems with a small number of components 0. 

On the other hand, colloidal suspensions are rather commonly polydisperse. Polydisper- 
sity means that mesoscopic suspended particles of a same chemical species are not necessarily 
identical, but some of their properties (size, charge, etc.) may exhibit a discrete or continu- 
ous distribution of values. Even when all macroparticles belong to a unique chemical species, 
a polydisperse fluid must therefore be treated as a multicomponent mixture, with very large 
p values - of order 10 1 4- 10 3 or more (discrete polydispersity) - or with p — > oo (continuous 
polydispersity). 

The above-mentioned shortcomings of the SHS1-PY solution offer a strong motivation 
for investigating an alternative sticky hard sphere model, proposed by Brey et al. [f| and 
Mier-y-Teran et al. ||, and hereafter referred to as SHS2 model. The adhesive part of its 
potential is defined starting from an attractive Yukawa tail, which, in the sticky limit, has 
both amplitude and inverse range tending to infinity, with their ratio remaining constant. In 
this case, the OZ equations are analytically solvable within the mean spherical approximation 
(MSA) §,[7j. Although the SHS2-MSA solution is simpler than the SHS1-PY one, it has 
received much less attention, especially in the multi-component case. 
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In two previous papers |8||| we investigated structural properties of polydisperse fluids 
using a version of the SHS2 model in which the coupling (stickiness) parameters, which 
define the strength of the Yukawa attraction, are factorizable. This choice is the simplest 
one 



myiTf . In fact, a slight different version of the SHS2 potential (with non-factorizable 
coefficients), proposed by Tutschka and Kahl ]T^-0 leads to more complicated analytical 
results, without any increase of the physical insight. 

Once that the structural properties of a model are known, the next natural step is to 
study the corresponding thermodynamics. Unfortunately, neither for SHS1 |§ nor for SHS2 
(see below), this is a simple task. 

In this short contribution we focus in particular on the compressibility equation of state 
(EOS) of the SHS2 multi-component model, since the compressibility route represents the 
simplest method for obtaining the pressure, given an analytical expression for Baxter's factor 
correlation function ?y(r). Under rather general conditions, we show that no compressibility 
EOS can exist for the SHS2 multi-component model, within the MSA. We argue how this 
inconsistency stems from the MSA closure, and show a possible way out to overcome this 
difficulty, by using a simple illustrative example. 



After this work has been completed, we became aware of Ginoza's recent analysis JT5 
where the author discusses a factorizable model essentially identical to the one considered 
in our previous work ||, but fails to recognize the MSA inconsistency, and thus reports an 
incorrect expression for the compressibility EOS. 



Our findings agree with those by Tutschka and Kahl [13,14], who observed the same 



inconsistency within their particular version of the SHS2 model. 

II. THE SHS2 MODEL 

The starting point of the SHS2 model is a fluid with particles interacting via a HS 
repulsion plus a particular Yukawa (HSY) attraction, i.e. 
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< SY (r) = { ' ' v ' ■ (1) 



< r < Gij = (o-j + (Tj)/2 



-zAijt z{j ai ^/r, r>aij. 

Here, cxj denotes the HS diameter of species % (whose number density is p^, z is the inverse 
range of the Yukawa tail, all Aij = Aji stickiness parameters are > 0, and the well depth, 
By = zAij/aij, depends on z linearly. The OZ equations for HSY mixtures have been solved 
analytically fi~6f , for any finite z, within the MSA closure, i.e., Cjj (r) = —f3uij(r) for r > 



(cij (r) is the direct correlation function, and (3 = (/cbT)^ 1 ). The MSA solution for SHS2 
can thus be obtained by taking the sticky limit, z — ► oo, of the solution gS SY_MSA (r) for the 
HSY fluid. The result is 

^ \di{r 2 - a?) + h(r- atj) + K ij7 L {j <r< a {j 

Qij[ r > — \ [*) 

0, elsewhere 

1 ggggj _ Xj 3^2g| X^ 

fli _ A + A 2 A ' 1 ~ 2A 2 + 2A ' 1 j 

7r 

Xj — ^ ' Pm^ra M-imi ^-im 12Xj m , (4) 

u m 

where Lij = (cr^ — cr,)/2, £ n = (vr/6) J2i Pi^n an d A = 1 — £ 3 . The coefficients = qij{cr^), 
given by 

7V-SHS2-MSA _ jjj /r\ 

are density-independent, and have dimensions of [length] 2 . 

Tutschka and Kahl's version of the SHS2 model [p!2|-[T4]1 hinges upon non-factorizable pa- 
rameters, t^shs2-msa _ j^ffi with 7^- obeying a Berthelot-type rule, i.e., 7^ = (ju^jj) 1 ^ 2 ■ 
Our version with factorized coefficients assumes that A^ = EoGiGj, where Eq is an energy 
and G m has dimensions of length. Thus 

1^ = , (6) 
where Y m = 7o G m , 7^ = = — L- , (7) 



with T* being a reduced temperature (as in Baxter's model, the factor 1/12 is introduced 
for later convenience. In our previous papers |8],|9| it was absent, and the correspondence 
between the two reduced temperatures is: T* = T* ld /12). 

Irrespectively of the choice for the coefficients Ay, the SHS1 and SHS2 potentials are 
different and should not be confused even in the sticky limit. An intuitive way of under- 
standing this point is to notice that, as the well width goes to zero, the area of the square 
well in Baxter's SHS1 model vanishes, whereas the area under the Yukawa tail in SHS2 
remains finite ||. This difference becomes important when evaluating thermodynamics. In 
fact, it can be shown that the virial pressure depends not only on the qij(r) resulting after 
the sticky limit, but also on the functional form of the tail as well as on the way the sticky 



limit is taken ||17|| . Furthermore, the SHS1 model is analytically solvable within the PY 
closure, but not within the MSA one, whereas the opposite is true for SHS2. 

At the level of approximate solution for ?y(r), the difference between the SHS2-MSA 
and SHS1-PY expressions lies only in the coefficients Ay, which in the latter case read 

7>-SHSl-PY 1 „.PY(„ \ 2 l x 2 /n\ 

K ij - TI7~: !Ji -i wfaj = 12 ( 8 ) 
where the dimensionless positive parameter ry, that appears in the SHS1 potential, is related 
to both the temperature and the stickiness between particles of species i and j, while yy(ay) 
is the contact value of the cavity function. Note that A^ HS1_PY is non-factorizable and 
density- dependent, since ?/y(cry) depends on the densities of all components in the mixture. 
This difference, albeit seemingly harmless, has far reaching consequences, as it will be shown 
in the following. 



III. COMPRESSIBILITY EQUATION OF STATE 

Once q%j{r) is known, one can calculate derivatives of the compressibility (c) pressure by 
means of two general relations, obtained from fluctuation theory in the grand-canonical en- 
semble and from Baxter's factorization of the three-dimensional Fourier transform of cy (r), 
i.e., 
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'df3P s 
, 9 Pi , 



= a,i - 27r^2p m a m q mi (0), (9) 

T,p k rn 



where qij (k) is the one-dimensional Fourier transform of q^ (r), while xt — P^bTKt de- 
notes the isothermal susceptibility [Kt being the isothermal compressibility), and dj = 
1 — 27T J2i PiQu(fy- Note that the pressure is a function of (T, p 1 , . . . , p p ) in Eq. @, and of 
(T, p, x) in Eq. fllCf), where x = (xi, . . . , x p _i) represents the composition in terms of molar 
fractions. 

We can apply the previous relations to SHS models for colloidal fluids, after observing 
that Eq. (|2]) yields q m i(0) = 12 _1 (a m af + 3A _1 cr m of — M mi Oi). Inserting this term and the 
a,i given by Eq. (^) into Eqs. (|9|)- ([TQ|) , one gets an expression for [d(3P '/ 'dpi) T pk and 

/^P\ fd(3P\ US ~ PYc 2(X) 66 (gjQ (X 2 ) 

I 9 P ; Tx ^ dp ; Tx A 2 A3 + A 2 ' 1 ) 

where we have introduced compositional averages, such as (/) = J2m x mfm and (fg) = 
J2m x mfm9m, and the derivative with superscript HS-PYc refers to the compressibility pres- 
sure of the corresponding HS mixture, evaluated within the PY approximation. Note that 
the two results hold for both the SHS2-MSA and SHS1-PY solutions, if the corresponding 
Kij's are used. 

In order to obtain the compressibility pressure, it might now be spontaneous to per- 
form immediately the integration of Eq. (|TTD with respect to p, at fixed composition x, 
i.e. PP(T, p, x) = JoX^dp ■ When applied to the first term on the r.h.s. of Eq. (PD, 
(dpP/dp)^ s ~ PYc , this procedure does indeed lead to the known PY C - EOS for HS mixtures. 
For the SHS models, the integration requires the knowledge of the dependence (if any) of 
the coefficients on density. 

In the SHS1-PY case, the aforesaid calculation is nevertheless practically impossible, 
due to the lack of an explicit expression for Ay(p). This difficulty has been by-passed by 
exploiting a further result by Baxter, which directly provides the PY pressure itself rather 
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than the inverse susceptibility @J3[]. On the other hand, in the SHS2-MSA case the density- 
independence of xfj 2 ~ might suggest that the integration of Xt 1 1S a straightforward 
operation. Unfortunately, this is not so. 

In fact, to ensure that the whole procedure is correct, one must first test whether the 
differential J2i (dpP/dpi) dpi, constructed with the partial derivatives given by Eq. (|9|), is 
exact, since P must be a state function. For this to occur, it is necessary that 

(12) 



d 


(d(3P\ 


d 


(d(3P\ 


d Pj 


[dp J 


dpi 


[ d Pj ) 



for any pair i and j. This symmetry condition, obeyed by any exact theory, may not be met 
when using an approximate closure. In this case, no compressibility EOS can exist within 
the considered approximate theory. Now, we show that the equality (|1^) of the mixed 
second-order partial derivatives is not necessarily satisfied for SHS mixtures. 

Let us assume, rather generically, that a certain SHS model with an appropriate closure 
has a solution %(r) given by Eqs. with density-independent, symmetric, coefficients 

Mij = VIKij = Mji. Inserting this solution into Eq. (§) and taking the derivative of d/3P/ dpi 
with respect to pj yields a rather lengthy expression for d 2 f3P/ dpjdpi. Upon discarding all 
terms which are apparently symmetric with respect to an exchange of indices i and j, we 
are left with the following sum 



S 1 (ij) = —GiMijXi + ^afajXj 

IT 



1 



of (X, - X]) + SafcXj + afa, (xf - - £ Pl X { M 3l 



(13) 



where xf^ = (tt/6) J2i PiMn. A compressibility EOS can exist, within an approximate theory 
of the considered kind, only if the relevant coefficients are such that Si(ji) = S\(ij). 



Tutschka and Kahl's choice |l2| , p^| does not meet this requirement. On the other hand, 
it is easy to verify that even any choice with factorized coefficients, = M GiGj, fails 
to satisfy the necessary condition (Mo is a density-independent factor, which in our model 
defined by Eqs. (§)-(0) coincides with 1/T*). S\(ij) cannot be symmetric for a generic choice 
of GmS, but the same occurs if we assume the power-law relationship between stickiness and 



size employed in Refs. [|||§, i.e. G m = / , where (a) denotes the average diameter 

and a > 0. In this case, the value a = 1 is however indicated as preferable, since it represents 
the only way of making the first term, — A~ l aiMijXj, symmetric. 

As a first result, one can thus state that no compressibility EOS can exist within the 
SHS2-MSA theory with coefficients given either by any factorization rule My = M^GiGj, or 
by Tutschka and Kahl's unfactorized choice |12|,[14[]. We note that in Ginoza's recent analysis 
P~5|1 , such a crucial feature is failed to be recognized, and the corresponding expression for 



the compressibility EOS is therefore incorrect. 

It is clear that the violation of the relation ( |12^ ) can be traced back to the inadequacy of 
the MSA closure, which is responsible for the density-independence of the coefficients 
(or Kij). More generally, one might suspect that no compressibility EOS can exist for any 
solution qij (r) given by Eqs. (0)-© as long as the M^-'s are density-independent. 

Let us now assume that the M^'s depend on the densities (pi, . . . , p p ). In this case, 
d 2 (3P / dpjdpi contains further contributions stemming from derivatives. Again discarding 
those that are clearly symmetric, one gets the following terms 

l r / J( 2 \dM u 

1 7r ^ ( „ , 6£ 2 q „ o \ dM ml 



l,m \ Jr. 

1 7T „ 8 (MuM m i) 



A c / j rlr id ^ id r\ 

6 U 8 pi 

2 3^\ 2 v- 8{M ml M nl ) 

l,m,n <3 

which must be added to Si(ij), to form a new symmetry condition for S(ij) = Si(ij)+S2{ij)- 
We emphasize the fact that the density- dependent coefficients of the SHS1-PY solution, 
as given by Eq. (||), should satisfy the requirement S(ji) = S(ij), but a direct test of 
this feature proves to be a highly non-trivial task. However, this observation prompts 
the suggestion that a simple density-dependent modification of the best SHS2-MSA choice, 
Mij = Mo<Ji<jj, along the lines of the SHS1-PY solution, might lead to a fulfilment of the 
symmetry condition. Indeed, we find that the choice 
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Mij = Mo^Wj , with M = (J^j , (15) 

(where is an arbitrary function vanishing as T* — > oo) leads to a compressibility pressure 
P satisfying the condition (|i2|). 

Next, we discuss a possible physical origin of the factor 1/A appearing in the above 
solution. First, we note that in Kf^ sl ~ PY of Eq. @ one could relate afj to a measure of 
adhesive interaction surface (since the centre of a particle i which moves around a particle 
j, but remains in contact with it, must lie on a spherical surface with radius cry), while the 
density-dependent factor yij(o~ij) represents the probability of finding a particle of species i 
touching any given particle of species j. On the other hand, in the corresponding modified 
expression ([15]), one could imagine a^aj as representing the area of an interaction spherical 
surface with radius (aiO-j) 1 ^ 2 , while 1/A may be reckoned as a crude approximation to Vij(o~ij) 
(1/A is indeed the simplest term, independent of the species indices i and j, appearing in 
the expression of cavity functions at contact). 

There are (at least) two ways to relate the result ( |ToD to some possible closure, which 
should represent an improvement over the MSA one. First, one might consider a generalized 
mean spherical approximation (GMSA) instead of the MSA. This amounts to replace 
the MSA, Cij(r) = —f3u^ SY (r) = z^A^e'^^ ^ jr for r > a^, with the Yukawa closure 
Cij (r) = zKije~ z( - r ~ ai ^ jr for r > cr^-, with parameters Kij not given by Eq. ([5]), but density- 
dependent and initially undetermined. The GMSA has often been used in the past, and its 
unknown coefficients have usually been determined by employing some thermodynamic 
consistency condition. In the present case, one could regard the symmetry condition (p~2|) 
as an alternative condition for the GMSA closure. Within this conceptual framework, Eq. 
(|T5|) is the simplest solution. 

Alternatively, we recall that the PY approximation may be written in the form Cy- (r) = 
fij ( r ) Uij ( r ) i while the MSA may be derived from it under the approximation f\j (r) ~ 
—(3uij(r), with (r) ~ 1, for r > cr^-. Hence one could define a new density-dependent 
closure for SHS potentials, assuming c^ (r) = —yij{cij) (3uij{r) for r > a^. The use of 



u 



HSY, 
ij I 



r) with factorized coefficients, along with the rough approximation yij{<Jij) — 1/A, 



would then lead to the M^- coefficients expressed by Eq. (0), with (1/T*) = 1/T*. 

Finally, it is worth reporting that the solution corresponding to Eq. ( p~5D yields the 
following compressibility EOS 

^ CV =I^ + ( 3 - M o)e 1 ( I ^) 2 +^(3-Mo) 2 e 2 ( T ^) 3 , (16) 

where V = (vr/6) (<r 3 ) denotes the average volume of a particle, rj = £3 = pV is the packing 
fraction, while the dimensionless parameters e\ = (cr) (a 2 ) / (cr 3 ), e 2 = (c 2 ) 3 / (<J 3 ) 2 depend 
on the molar composition x and reduce to e\ = e 2 = 1 for one-components fluids. Note that 
M = (p (1/T*) may or may not coincide with 1/T*, but it must vanish as T* — > 00, so that 
the PYc - EOS for HS mixtures is recovered. 

The simplicity of Eq. flT6|) is due to the fact that now Xj = M ^2C"«/A, and thus we get 

a, = ^ + (3-M )^, & l = -(3-M )|g|, (17) 

which differ from their HS counterparts only in having the HS factor 3 replaced with the 
temperature-dependent coefficient 3 — M . 

In spite of its plainess, Eq. (|T6"D represents an analytical and consistent compressibility 



EOS for a multicomponent system of particles with both repulsive and attractive interac- 
tions. Furthermore, it bears interesting connections with the SHS1-PY model which will be 
discussed elsewhere. 



IV. CONCLUSIONS 

The PY compressibility EOS of Baxter's SHS1 model is practically inapplicable to fluids 
with a large number of components. This difficulty urges to search for either alternative 
closures or different models. 

In this paper we have investigated the SHS2 model, which is analytically solvable within 
the MSA, and whose thermodynamics is still rather unexplored, especially in the multi- 
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component case. In particular we have focused on the compressibility EOS, since it is the 
simplest route to get the pressure from q^ (r). 

Mixtures require a more careful analysis than pure fluids, and one should always consider 
the possibility that an approximate theory suffers from some thermodynamic inconsistencies 
not present in the one-component counterpart. In particular, in order to ensure the existence 
of the compressibility EOS, we have pointed out the necessity of an explicit check of a basic 
thermodynamic consistency requirement hinging on the equality of the mixed second-order 
partial derivatives of P with respect to the densities, which , in turn, can be expressed in 
terms of (r) . This crucial feature is often overlooked in the literature, perhaps because it 
would be automatically fulfilled in an exact theory. In an approximate theory, on the other 
hand, this is not the case, and the reciprocity condition ( |P^ ) is one of the sum-rules that 
must be explicitly checked for. 

Two are the main results of the present work. A first finding is that, although an MSA 
compressibility EOS is known for pure SHS2 fluids |J, its extension to mixture is not pos- 
sible, since the SHS2-MSA solution violates the aforesaid symmetry condition. As a direct 
consequence, the compressibility EOS recently reported by Ginoza |THj is flawed. We have 



discussed how this inconsistency occurs for any choice of factorized stickiness coefficients, as 
well as for Tutschka and Kahl's unfactorized ones [ 12jl4| , and argued that it has its origin in 



the deficiency of the MSA closure. The reason appears to be the density-independence of the 
MSA coefficients = qij(cr^). Second, we have presented an illustrative example, where 
the inclusion of a simple and plausible density- dependence in the matrix Ky, produces a 
fulfilment of the required condition, and generates a possible compressibility EOS. 

This result seems to suggest that no compressibility EOS can exist for any closure lead- 
ing to a solution q^ (r) of the form given by Eqs. (@)-(f§) with density- independent 
coefficients. This also prompts the necessity of abandoning the MSA and resort to density- 
dependent closures. We have attempted to do this by interpreting our result for the com- 
pressibility EOS, as originating from a different, more sophisticated, approximate theory. 

Although it is clear that our result cannot be considered as the correct final solution to 
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our problem, it nevertheless represents, not a simple academic exercise, but a useful step 
towards a satisfactory EOS for SHS mixtures. 
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